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ABSTRACT 

We review the status of the neutrino oscillation physics (as of June 2003), with 
a particular emphasis on the present knowledge of the neutrino mass-mixing pa- 
rameters in a three generation approach. We consider first the — ► v T flavor 
transitions of atmospheric neutrinos. It is found that standard oscillations provide 
the best description of the SK+K2K data, and that the associated mass-mixing 
parameters are determined at ±ler (and Adf = 1) as: Am 2 = (2.6 ± 0.4) x 10 -3 
eV 2 and sin 2 29 = l.OOloios- Such indications, presently dominated by SK, could be 
strengthened by further K2K data. Then we analyze the energy spectrum of reactor 
v events recently observed at KamLAND and combine them with solar and terres- 
trial v data. We find that the solution to the solar v problem at large mixing angle 
(LMA) is basically split into two sub-regions, that we denote as LMA-I and LMA- 
II. The LMA-I solution, characterized by lower values of the squared neutrino mass 
gap, is favored by the global data fit. Finally, we briefly illustrate how prospective 
data from the SNO and KamLAND can increase our confidence in the occurrence of 
standard matter effects in the Sun, which are starting to emerge from current data. 
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1 Introduction 



In its first phase of operation (years 1996-2001), the Super-Kamiokande (SK) exper- 
iment has provided, among other important results, compelling evidence for atmo- 
spheric disappearance PHI2]- This evidence, now firmly based on a high-statistics 
92 kton-year exposure [3] , has not only been corroborated by consistent indications 
in the MACRO [I] and Soudan 2 5J atmospheric neutrino experiments, but has 
also been independently checked by the first long-baseline KEK-to-Kamioka (K2K) 
accelerator experiment jHllZj, using SK as a target for v^s produced 250 km away 
with (E u ) ~ 1.3 GeV. Neutrino flavor oscillations, interpreted in terms of nonzero 
mass-mixing parameters (Am 2 , sin 2 29) in the — > u T channel, provide by far the 
best and most natural explanation for the observed disappearance PUI2]- 

In Section 2 we review the phenomenological status of the standard os- 
cillations in the — > v T channel, in the light of the latest SK atmospheric zenith 
distributions j3] and of the first spectral results from the K2K experiment [7j. 

On the solar neutrino front, the Chlorine [S], Gallium 01101 HI], Super- 
Kamiokande (SK) US ECU and Sudbury Neutrino Observatory (SNO) [H 113 EE] 
solar neutrino experiments have convincingly established that the deficit of the ob- 
served solar v e flux with respect to expectations [T7] implies new neutrino physics. 
In particular, the charged and neutral current (CC and NC) data from SNO have 
proven the occurrence of v e transitions into a different active state v a with a statis- 
tical significance greater than 5a [15J. 

Barring sterile neutrinos and nonstandard v interactions, such transitions 
can be naturally explained by the hypothesis of flavor oscillations i n the v e — + v a 
channel [v a being a linear combination of and v T ) driven by nonzero v squared 
mass difference and mixing angle parameters (5m 2 , # 12 ) JH]. The (v^, v T ) combi- 
nation orthogonal to v a is probed by atmospheric v oscillations [2 , with different 
parameters (Am 2 , #23) The third mixing angle 13 , needed to complete the 3x3 
mixing matrix, is constrained to be small by additional reactor results 122] , and 
can be set to zero to a good approximation for our purposes. 

The recent results from the Kamioka Liquid scintillator AntiNeutrino De- 
tector (KamLAND) [231 have provided a beautiful and crucial confirmation of the 
solar v e oscillation picture through a search for long-baseline oscillations of reactor 
!7 e 's. The observed of V e disappearance in KamLAND has confirmed the previ- 
ously favored solution in the (5m 2 , # 12 ) parameter space, often referred to as the 
large mixing angle (LMA) region [16J in the literature (see, e.g., [21] and references 
therein). Moreover, the KamLAND data have basically split this region into two 
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allowed subregions, which we will refer to as LMA-I and LMA-II, following Ref. |25j . 

In Sections 3 we analyze the first KamLAND spectral data [23] and com- 
bine them with current solar neutrino data (HI El E3 El Ell El El E], assuming 
two-flavor oscillations of active neutrinos [21], in order to determine the surviving 
sub-regions of the LMA solution. In the analysis we include the CHOOZ reactor 
data [21]. 

Finally, in Section 4 we briefly illustrate how emerging indications of solar 
matter effects can be corroborated in the LMA parameter region. In particular, we 
show that the amplitude of matter effects (introduced as a free parameter omsw) can 
be significantly constrained by using prospective data from SNO and KamLAND. 

2 "Atmospheric" neutrinos 

A careful analysis of the SK and K2K data sets used in the following can be found in 
j2H]- Concerning SK atmospheric neutrino data (92 kton-year [3]), we use the usual 
zenith angle (9 Z ) distributions of leptons: sub-GeV e-like and /i-like events, divided 
in 10+10 bins; multi-GeV e-like and /i-like events, divided in 10+10 bins; upward 
stopping and through-going /i events, divided in 5+10 bins. The calculation of the 
theoretical events rates i?^ heo in each of the 55 bins is done as in [2*71 l2*Sl HQ . The 
SK statistical analysis is considerably improved with respect to [23 [2H]- Now the 
set of systematic errors has been enlarged to 11 entries, leading to a more complex 
structure of correlated errors affecting the i?Jf ec "s, as emphasized in |3*U] . 

Concerning the K2K data, we use the absolute spectrum of muon events 
in terms of the reconstructed neutrino energy E [7J, which provides a total of 29 
events (here divided in 6 bins). In this sample, the parent neutrino interactions 
are dominantly quasi-elastic (QE), and the reconstructed energy E is thus closely 
correlated with the true neutrino energy E v . 

Let us now discuss the updated bounds on the parameters (Am 2 , sin 2 29), 
governing the scenario of standard oscillations (here 9 = #23)- 

Fig. 1 shows the joint bounds on the (Am 2 , sin 2 29) parameters from our 
analysis of SK, K2K, and SK+K2K data. The bounds in the left panel are very 
close to the official SK ones, as presented in [3]. The bounds in the middle panel are 
instead slightly weaker than the official K2K ones [7], especially in terms of sin 2 29. 
In particular, we do not find a lower bound on sin 2 29 at 99% C.L. (for N^-p = 2). 
The reason is that we cannot use the additional (dominantly) non-QE event sample 
of K2K (27 events), which would help to constrain the overall rate normalization 
and thus sin 2 29. This fact might also explain why we find the K2K best fit at 
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Figure 1: Standard oscillations in the — > v T channel: bounds on the parameters 
(Am 2 , sin 2 26) from SK atmospheric data (left panel), K2K spectral data (middle 
panel), and their combination (right panel). 



sin 2 26 = 0.82 rather than at 1.00 as in jjj. By comparing the left and right panels 
of Fig. 1, the main effect of K2K appears to be the strengthening of the upper 
bound on Am 2 , consistently with the trend of the first K2K data (rate only no 
spectrum) 29J. The main reason is that, for Am 2 ~ (4-6) x 10 -3 eV 2 , the first 
oscillation minimum would be located at — or just above — the K2K energy spectrum 
peak, implying a strong local and overall suppression of the expected events. 

Fig. 2 shows on the left the SK and SK+K2K bounds on Am 2 , when the 
sin 2 26 parameter is projected (minimized) away. The linear scale in Am 2 makes the 
K2K impact on the upper limit more evident. Notice that, up to ~3cr, the global 
(SK+K2K) x 2 function is approximately parabolic in the linear variable Am 2 , so 
that one can define a one- standard-deviation error for this parameter. This feature 
was already argued on the basis of a graphical reduction of the official SK and K2K 
likelihood functions [21], and is here confirmed through a full analysis. By keeping 
only the first significant figure in the error estimate, a parabolic fit provides the ±1<t 
range 



Am 2 



v. 



(2.6 ± 0.4) x 10" 3 eV 2 . 

The bounds on sin 2 26 are instead entirely dominated by SK. This is shown 
on the right of Fig. 3, where the A\ 2 function in terms of sin 2 26 is reported, for 
Am 2 projected (minimized) away in the SK fit. Here the addition of K2K data 
would insignificantly change the bounds (not shown), which thus hold for both the 
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Figure 2: Standard oscillations in the — > v T channel. On the left: bounds on Am 2 
for unconstrained sin 2 26 from SK (dashed curve) and SK+K2K (solid curve). On 
the right: bounds on sin 2 26 for unconstrained Am 2 from SK data. The inclusion of 
K2K data induces here negligible changes (not shown). 

SK and the SK+K2K fit. Also in this case, the nearly parabolic behavior of Ax 2 
allows to properly define a la range, 

sin 2 26 = 1.00±g;gg , (2) 

with the lower Na error scaling linearly with N (up to iV ~ 3). Equations (|T|) 
and (J2J concisely review the current fit to the standard oscillation parameters, as 
anticipated in the Introduction. 

Fig. 3 shows the comparison between observations and best-fit predictions 
for the SK zenith distributions. In particular, the comparison between solid and 
dashed histograms shows that the systematic shifts are often comparable in size to 
the statistical errors, implying that just increasing the SK atmospheric v statistics 
will hardly bring decisive new information on the standard oscillation scenario. In 
the SG and MG samples, the fit clearly exploits the systematic uncertainties to 
increase the e-like event normalization, especially in the upward direction, so as to 
reduce the "electron excess" possibly indicated by SK data. 

Concerning //-like events in the SG and MG samples, the fit shows an 
opposite tendency to slightly decrease the normalization of (especially down-going) 
events. The tendency appears to be reversed in the high-energy UT sample. Taken 
together, these opposite shifts of e-like and //-like expectations in the SG and MG 
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Figure 3: SK experimental zenith distributions (ztl<7 sta t,), compared with the corre- 
sponding theoretical ones at the global (SK+K2K) best-fit point. All distributions are 
normalized to the unoscillated predictions in each bin. For the theoretical event rates, 
we show both the central values R^ eo (dashed histograms) and the "shifted" values 
R n (solid histograms), which embed the effect of systematic pulls. The difference 
between R n and i?^ heo shows how much (and in which direction) the correlated 
systematic errors tend to stretch the predictions in order to match the data. 



samples seem to suggest some systematic deviation from the predicted fi/e flavor 
ratio which, although not statistically alarming, should be kept in mind: deviations 
of similar size might have their origin in neutrino physics beyond 2u oscillations. 
Unfortunately, since such effects are typically not larger than the systematic shifts 
in Fig. 3, they are likely to remain hidden in higher-statistics SK data. 



3 Impact of KamLAND on solar neutrinos (a 2u analysis) 



The KamLAND recent observation of V e disappearance j23j confirms the current 
interpretation of solar neutrino data [T3J EH "2H ED] in terms of v e — > v^ T oscilla- 
tions induced by neutrino mass and mixing [Tj"J [""J"] , and restricts the corresponding 
parameter space (5m 2 , 812) within the so-called large mixing angle (LMA) region. 
In this region, globally favored by solar neutrino data, matter effects |31| 13*2] in 
adiabatic regime [3*3*] I3T] are expected to dominate the dynamics of flavor transi- 
tions in the Sun (see, e.g., [S]). The KamLAND spectral data appear to exclude 
some significant portions of the LMA solution [23], where the predicted spectrum 
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Figure 4: Two-flavor active neutrino oscillations: Global analysis of solar and 
CHOOZ neutrino data in the (5m 2 , sin 2 # 12 ) parameter space, restricted to the LMA 
region. The best fit is indicated by a black dot. 

distortions [21] would be in conflict with observations [2Hj . 

In the 2z/ case, we find that the inclusion of the KamLAND spectrum 
basically splits the LMA solution into two sub-regions at "small" and "large" 5m 2 , 
which we call LMA-I and LMA-II, respectively (the LMA-I solution being preferred 
by the data) |35| . Such regions are only slightly modified in the presence of 3u 
mixing, namely, for nonzero values of the mixing angle #13. We also present updated 
bounds (as of June 2003) in the 3z/ parameter space (5m 2 , #12, #13). 

In our KamLAND data analysis |35j . we use the absolute spectrum of 
events reported in [23], taken above a background-safe analysis threshold of 2.6 MeV 
in visible energy E. The events below such threshold might contain a significant 
component of geological I7 e 's |36| . whose large normalization uncertainties are poorly 
constrained at present by the KamLAND data themselves (221- Above 2.6 MeV, a 
total of 54 events is found, against 86.8 events expected from reactors [23 ■ Finally, 
the observed energy spectrum of events is analyzed as in [21] , with the improvements 
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Figure 5: Two-flavor active neutrino oscillations: Analysis of the KamLAND energy 
spectrum data above 2.6 MeV in the (5m 2 , sin 2 # 12 ) parameter space. A "tower" of 
octant-symmetric regions is allowed at different values of 5m 2 . The symmetric best 
fits are indicated by black dots. The left dot is remarkably close to the solar best fit 
in Fig. 4- 

reported in [35j. 

An updated 2u analysis of current solar+CHOOZ neutrino data, [21], is 
presented here. The fit includes 81 solar neutrino observables [301121] an d 14 CHOOZ 
spectrum bins |21| I24|. for a total of 95 data points. The A% 2 expansion around the 
minimum, relevant for the estimation of the oscillation parameters (5m 2 , sin 2 6 I 12 ), is 
shown in Fig. 4, where we have restricted the 5m 2 range to the only three decades 
relevant for the LMA solution and for the following KamLAND analysis. 

In Fig. 5 we report the 2z/ analysis of KamLAND [35J: there appears 
to be a "tower" of solutions which tend to merge and become indistinguishable 
for increasing 5m 2 ; the lower three ones are, however, rather well separated at 
90% C.L. Notice that our allowed regions are slightly larger (i.e., less constraining) 
than those in the rate+shape analysis of [23]- One of the two octant-symmetric 
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Figure 6: Two-flavor active neutrino oscillations: Global analysis of solar, CHOOZ, 
and KamLAND neutrino data in the (5m 2 , sin 2 12 ) parameter space. With respect 
to Fig. 4, the LMA region is significantly restricted, and appears to be split into two 
sub-regions (LMA-I and LMA-II), well-separated at 99% C.L. 



best fits points in Fig. 5 (black dots) is remarkably close to the best fit in Fig. 4. 
The difference in location with respect to the KamLAND official best-fit point at 
(Sm 2 /eV 2 , sin 2 9 12 ) = (6.9 x 10~ 5 , 0.5) [23] is not statistically significant, amounting 
to a variation A% 2 <C 1. 

The combination of the solar+CHOOZ results in Fig. 4 with the Kam- 
LAND results in Fig. 5 gives the global 2u results shown in Fig. 6. Two rather 
distinct solutions, that we label LMA-I and LMA-II, are seen to emerge. The LMA- 
I solution is clearly preferred by the data, being close to the best fit points of both 
solar+CHOOZ and KamLAND data. The LMA-II solution is located at a 5m 2 
value about twice as large as for the LMA-I, but is separated from the latter by a 
modest A% 2 = 5.4 difference (dominated by solar neutrino data). Indeed, if we con- 
servatively demand a 99.73% C.L. for the allowed regions, the LMA-I and LMA-II 
solutions appear to merge (and extend towards Sm 2 ~ 3 x 10~ 4 eV 2 ) in a single 
broad solution. In any case, at any chosen C.L., the allowed regions of Fig. 6 are 
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significantly smaller than those in Fig. 4. Therefore, with just 54 initial events, the 
KamLAND experiment is not only able to select the LMA region as the solution to 
the solar neutrino problem, but can also significantly restrict the corresponding os- 
cillation parameter space. With several hundred events expected in the forthcoming 
years, there are thus very good prospects to refine the parameter estimate [21] . 

4 Indications of matter effects in the Sun 

Within the LMA region, solar neutrino oscillations are governed not only by the kine- 
matical mass-mixing parameters (5m 2 , 812), but should also be significantly affected 
by the interaction energy difference (V = V e — V a ) between z/ e 's and u a 's propagating 
in the solar (and possibly Earth) background matter |HUIS21, through the so-called 
Mikheev-Smirnov-Wolfenstein (MSW) mechanism [3T] in adiabatic regime (SHI- Al- 
though Earth matter effects (i.e., day-night variations of solar event rates) remain 
elusive, solar matter effects seem to emerge, at least indirectly, from the combina- 
tion of the available data (and especially from SNO), through a preference for an 
average oscillation probability smaller than 1/2 at energies of a few MeV. A phe- 
nomenological approach to the problem has been recently presented in [33] , where a 
free parameter czmsw is introduced, called to modulate the overall amplitude of the 
interaction energy difference V in the dynamical term 7id y n of the Hamiltonian, 

V -> a M sw • V . (3) 

By treating czmsw as a continuous parameter, one can try to constrain its allowed 
range through global data analyses: A preference for omsw ~ 0(1) would then 
provide an indirect indication for the occurrence of matter effects with standard 
size, as opposed to the case of pure "vacuum" oscillations (Hd yn — 0). 

We have verified that the current solar neutrino data, by themselves, place 
only very loose and uninteresting limits on omswj as f ar as the mass-mixing oscilla- 
tion parameters are left unconstrained. In fact, since the oscillation physics depends 
mostly on the ratio V/k, where k = 5m 2 /2E is the neutrino oscillation wavenumber, 
a variation of the kind V —>■ czmsw^ is largely absorbed by a similar rescaling of 
k (i.e., of 8m 2 ). In order to break this degeneracy, we need to include explicitly 
an experiment which is highly sensitive to 5m 2 and basically insensitive to matter 
effects, such as KamLAND. 

A SSM-independent preference for P ee < 1/2 has been provided first by 
the combination of SNO CC and SK data jHj and then by SNO data alone through 
the CC/NC double ratio jTH], but not yet with a significance high enough to rule out 
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P ee = 1/2 [3U]. Let us consider, in particular, the latest SNO constraints in the plane 
($ e , charted by the solar u e and v^ T fluxes, as shown in Fig. 3 of the original 
SNO paper ^5]. In such a figure, although the SNO best-fit point clearly prefers 
P ee ~ 1/3 (corresponding to $ MT ~ 2$ e ), the 95% C.L. ellipse is still compatible 
with P ee ~ 1/2 (namely, $ MT ~ $ e ). However, future SNO NC and CC data can 
considerably improve the constraints on P ee , by reducing both the statistical and 
the systematic error on the CC/NC ratio [3*7] . 

In conclusion, although the combination of all current solar neutrino data 
suggests a pattern of P ee compatible with the LMA energy profile and indicates an 
overall preference for the first octant of 9i 2 [16| . the emerging indications in favor of 
solar matter effects from this data set are not strongly compelling yet. 

Until now we have illustrated how a single datum (the SNO CC/NC double 
ratio) can discriminate the case of standard matter effects (omsw = 1) from the 
case of zeroed matter effects (omsw — 0) in the LMA parameter region. By using 
further experimental information from KamLAND, one could try to test [31] whether 
the "solar + KamLAND" combination of data can constrain matter effects in the 
Sun to have the right size [czmsw ~ ^(1)]- I n this kind of analyses, KamLAND 
basically fixes the oscillation parameters (5m 2 , sin 2 #12), and thus the kinetic part 
of the Hamiltonian, 7ikin- The role of solar neutrino data is then to check that the 
overall amplitude cimsw of the interaction energy difference V in the dynamical term 
Hdyn is consistent with the standard electroweak model (cxmsw — !)• 

We have thus performed global analyses including both current solar neu- 
trino data and current (or prospective) KamLAND data, with (5m 2 , sin 2 #12, cimsw) 
unconstrained. In particular, the analysis of current KamLAND data is based on 
the binned energy spectrum of reactor neutrino events observed above 2.6 MeV (54 
events) j2H|- Prospective KamLAND spectral data have instead been generated, 
with the same energy threshold and binning, by assuming either the LMA-I best-fit 
point (5m 2 = 7.3 x 10~ 5 eV 2 and sin 2 12 = 0.315) or the LMA-II best-fit point 
(5m 2 = 15.4 x 10~ 5 eV 2 and sin 2 9 12 = 0.300) [22], and increased statistics (5 x 54 
and 10 x 54 events). The CHOOZ reactor data [21] are also included. 

Figure 7 shows the results of such global fits, in terms of the function 
— X 2 ~ Xmm f° r variable omsw and unconstrained (i.e., minimized away) mass- 
mixing parameters. The na bounds on omsw are then given by Ax 2 = n 2 . Let us 
focus first on the solid curve, which refers to the fit with current KamLAND data. 
It appears that such curve can already place > 3cr upper and lower bounds on omsw- 
In particular, the hypothetical case of zeroed matter effects is already disfavored at 
~ 3.5cr, thus providing an indirect indication in favor of matter effects in the Sun. 
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Bounds on a MSW (Solar + CHOOZ + KamLAND) 




Figure 7: Bounds on omsw f or unconstrained (5m 2 , sin 2 6 , 12 ), including current solar 
and CHOOZ neutrino data, as well as current or prospective KamLAND data. The 
solid curve refers to the fit including current KamLAND spectrum data above 2.6 
MeV threshold (54 events), and shows that the hypothetical case of zeroed matter 
effects is already disfavored. The other curves refer to simulated KamLAND data, 
generated by assuming the LMA-I solution, and statistics increased by a factor of 
five (dotted curve) and of ten (dashed curve). 



The best-fit value of omsw is close to the standard prediction (omsw = 1)- However, 
the overall ±3<r range for aMsw, spanning about three orders of magnitude, is rather 
large. The width of this range can be understood by recalling the following facts: 
(1) the LMA range of 5m 2 constrained by solar neutrino data, which covers about 
one decade [3UJ 121]; can be shifted up or down by shifting omsw with respect to 
1, since the LMA oscillation physics depends on V/k oc aMsw/5m 2 ; (2) the range 
of Sm 2 constrained by current terrestrial data (including KamLAND+CHOOZ), 
which covers about two decades [22], is much less affected by omsw variations. 
As a consequence, by appropriately shifting a MSW , it is possible to overlap the 
reconstructed ranges of Sm 2 from solar and from reactor data over about 1 + 2 
decades. When the overlap sweeps through the degenerate 5m 2 intervals allowed 
by KamLAND alone [22], the fit is locally improved, leading to a "wavy" structure 
in the Ax 2 . In conclusion, although current solar+reactor data strongly disfavor 
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cfMSW — (zeroed matter effects) and provide a best fit close to omsw — 1 (standard 
matter effects), the presence of other local minima in the A% 2 function, as well as 
the broad 3er allowed range for OmsW; do not allow to claim a clear evidence of 
standard matter effects from current data. 

The broken curves in Fig. 7 refer to prospective KamLAND data, generated 
by assuming as true solution the LMA-I best-fit point. The energy threshold, the 
binning, and the systematic uncertainties are assumed to be the same as for the 
current KamLAND data. The dotted (dashed) curve refers to a number of reactor 
neutrino events five (ten) times larger than the current statistics. It can be seen 
that the global fit will progressively constrain omsw within one decade at ±3<r and, 
most importantly, will lead to a marked preference for cxmsw — 1) which is not yet 
evident in the present data. In conclusion, if the LMA-I solution is the true one, 
there are good prospect to test unambiguously the occurrence and size of standard 
matter effects in the Sun. 

5 Conclusions 

We have analyzed in detail the current SK atmospheric neutrino data and and 
the first K2K spectral data, in order to review the status of standard — > v T 
oscillations. We have then provided updated bounds for the standard oscillation 
parameters. In particular, the statistical analysis of the uncertainties reveals that 
K2K will lead further progress in this field, especially through higher-statistics tests 
of the low-energy spectrum bins. 

Going to solar neutrinos, the KamLAND experiment has clearly selected 
the LMA region as the solution to the solar neutrino problem, and has further 
reduced the (5m 2 , sin 2 # 12 ) parameter space for active neutrino oscillations. In the 
2v case, we find that the post-KamLAND LMA solution appears to be basically 
split into two sub-regions, LMA-I and LMA-II. The LMA-I solution, characterized 
by 5m 2 ~ 7 x 10 -5 eV 2 and sin 2 9 12 ~ 0.3, is preferred by the global fit. The LMA-II 
solution represents the second best fit, at about twice the value of 5m 2 . 

In the simplest picture, solar neutrino oscillations depend on the kinemat- 
ical parameters (5m 2 , sin 2 9 12) and on standard dynamical MSW effects in matter. 
These effects in current solar neutrino data are starting to emerge through an in- 
creasingly marked preference for P ee < 1/2, but still remain not clearly identified. In 
order to quantify statistically the occurrence of MSW effects, we have introduced a 
free parameter omsw modulating the amplitude of the v interaction energy difference 
in the neutrino evolution equation. By treating omsw as a continuous parameter, we 
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have then performed a global analysis including current solar, CHOOZ, and Kam- 
LAND data. The results are encouraging, since upper and lower bounds on omsw 
appear to emerge at the > 3a level. In particular, the case of "zeroed" matter 
effects is significantly disfavored. Moreover, the best-fit is tantalizingly close to the 
standard expectations for matter effects (cxmsw = !)■ 

The situation will improve through higher KamLAND statistics. In both 
the LMA-I and LMA-II cases, it appears possible to reduce the current uncertainty 
on a MSW by about two orders of magnitude. In conclusion, the selection of a single 
solution in the LMA oscillation parameter space appears to be crucial, before any 
definite conclusion can be made on the emerging indications of standard matter 
effects in the Sun. 
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